Phase behaviour of a hard sphere colloidal system in the presence 

of an external laser field 



Pinaki Chaudhuri, Ajay K. Sood, H. R. Krishnamurthy 
Department of Physics 
Indian Institute of Science 
Bangalore 560012, India. 



Abstract 



We study the phase behaviour of a two-dimensionally confined hard sphere 
colloidal system in the presence of a periodic light field of two interfering laser 
beams using Monte Carlo simulations. For a given packing fraction of the 
particles, the colloidal system undergoes a transition from a modulated liquid 
to a modulated crystal as the light intensity is increased, corresponding to 
laser induced freezing. For certain packing fractions of the colloidal particles, 
the system again becomes a modulated liquid as the field strength crosses a 
threshold value, showing a re-entrant behaviour. 
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In recent times, there has been a renewed interest in laser induced freezing (LIF), i.e, 
the freezing of colloidal particles in the presence of an external modulating (laser) field. 
LIF was first demonstrated by Choudhury et al [JTJ] . A sample of 1/im diameter polysterene 
particles was suspended in water placed between two optical flats. The highly charged 
particles, repelled by the walls, formed a two-dimensional monolayer structure. This two- 
dimensional system was subjected to a laser intensity pattern periodically modulated along 
one direction. It was found that when the wavevector of the stationary electric field is tuned 
to half the wave vector q Q where the liquid structure factor peaks, a triangular lattice with 
full two-dimensional symmetry results for light intensities above a threshold value. 

A simple theoretical analysis in terms of the Landau-Alexander-McTague theory 
showed a rich phase diagram with a first order freezing transition turning to second order 
at higher field strenghts via a tricritical point. In a density functional theory || analysis 
Chakrabarti et al |4[] obtained a similar cross over and in addition deduced a symmetry 
criteria for obtaining the tricritical point in the phase diagram. Monte Carlo simulational 
studies of 1.1/xm diameter colloidal particles interacting via a purely repulsive DLVO poten- 
tial by Chakrabarti et al || semed to confirm such a phase diagram, showing in addition an 
interesting reentrant transition to a liquid state for higher laser fields. Recently Wei et al 
M reported an experimental study of LIF and from the real space density plots concluded 
that there is indeed such a re-entrant liquid phase. However, more recent Monte-Carlo sim- 
ulations by Das et al |7||| of the same system as in || did not show any re-entrant phase. 
An even more recent theoretical work by Frey et al 0, argues (by extending the KTHNY 
picture of unbinding of dislocation pairs during melting in two dimensions to the case when 
an external field is present) that for highly-screened short-range interactions, there should 
exist a re-entrant transition to a liquid-state. 

In this paper, we study the phase behaviour of a system of hard sphere colloidal particles 
confined to a 2-d layer under the influence of a modulated laser field as a function of the light 
intensity. Using Monte Carlo simulations, we see that as the field intensity is increased, the 
system changes from a modulated liquid to a crystal, and finally back to a modulated liquid 
again, i.e, LIF in a hard-sphere colloid system does seem to show reentrant laser induced 
melting transition (LIM) unlike the 1.1/im DLVO colloidal system 

Hard sphere systems are interesting model systems in their own right [TIJ and should be 
a reasonable zeroth approximation for describing charge stabilised colloids which are highly 
screened. In such systems, where there is an exclusion of the interpenetration of the particles, 
the inter-particle potential is U(r) = oo(r < R) and U(r) = 0(r > R) , where R is the hard 
sphere radius. The free energy is determined entirely by entropy, which is dependent on 
the fraction of the total volume occupied by the spheres. At low volume fractions, the 
inter-particle collisions are rare. However when the packing fraction increases, the motion 
of the particles is restricted because of more and more collisions with neighbouring particles. 
As (f> is increased from the liquid phase, it becomes entropically favourable for the system to 
form a periodic crystal rather than a random fluid structure. Thus, there is a purely entropy 



driven liquid-to-solid transition which is interesting [ IT5| . The objective of this paper is to 
study the effect of the external field on this entropy driven phase transition. 

Our simulations have been done in an ensemble where both the rectangular box size 
and the number of particles are fixed. In order to change the the area fraction occupied 
by the hard spheres, the diameter of the colloidal particles is changed. Temperature is 
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not a relevant quantity when the external field is absent; but in the presence of the field, 
the unit of energy chosen is the thermal energy k^I '. We consider the external laser field 
to be modulated along the y-direction and model its effects in terms of a potential of the 
form V(r) = V cos(q y), where q a , equal to 2n/ 9 ^- ) is the smallest reciprocal lattice vector 
corresponding to the triangular lattice of lattice contant a. 

The particles are contained in a rectangular box where the ratio of the two sides L y 
and L x is : 2. This ensures that the simulation box is commensurate with a triangular 
lattice. L x (= L) is chosen as the square root of the number of particles. To avoid surface 
effects, we have chosen periodic boundary conditions using minimum image convention. 
The rectangular simulation box is replicated throughout space and in the course of the 
simulation, as a particle moves in the original box, so do the periodic images in exactly the 
same manner. In our similations, we have considered a system of 900 particles. The typical 
number of configurations over which equilibration has been done is of the order of 5 x 10 6 . 

The thermodynamic quantity of our interest is the translational order parameter. The 
translational order parameters are defined by : 



Pd k =< Y, ex P( i ^k-r l ) > (1) 

i 

In the above equation, is the real-space co-ordinate of the 2-th particle, q& is the k-th 
reciprocal lattice vector (rlv) belonging to the first shell of the rlvs. The averaging has been 
done over a canonical ensemble. 

The value of p^ k as defined in eqn (1) depends on the origin of co-ordinates. Therefore, in 
our simulations, we measure the translational order parameter as : 



Pq fc =< jr\ J2 ex P(iqk.ri)\ > (2) 

i 

The order parameter so defined is of order unity in the crystalline state, while in the liquid 
state it goes to zero as ^= for a system of N particles. 

The primitive translation vectors of the direct triangular lattice can be chosen as : 
ai = a(l,0),a 2 = a(-^r, 2 ) an< ^ those for the reciprocal triangular lattice as qi = 
0^/2 ( 2 ' 2O' ^ 2 = avI/2 ^' ^' ^ n our s i mu l a tions, we have chosen a to be unity. We denote 
by pi the order-parameters for the wavevectors ±q 2 parallel to the modulation wave-vector 
and pa denote the order parameter for the other four first-shell wave-vectors. 

In Fig.|l](a) we have plotted pd for the cases when there is no external field and when 
the field strength is infinite. In the case of zero field, entropy is the only mechanism that is 
forcing the density modes at the wave vectors to develop, and this happens at a packing 
fraction of about 0.7. For the case of infinte strength, the external potential is so large that 
the freedom of the particles to move in a y-direction is totally supressed and hence they can 
execute motion only in the x-direction. Now we find that pd becomes sizeable when the area 
fraction exceeds 0.68. This value of the area-fraction corresponds to the situation when there 
is vertical contact between the particles for the first time. For all area fractions less than 
0.69, the system consists of L uncoupled 1-dimensional chains of particles, with particles 
in the same line interacting through the hard-sphere repulsive potential. We know that 
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a disorder-to-order transition is not possible in a 1- dimensional system with short-ranged 
interactions. Hence there is no ordering before cj> = 0.69. Development of the density mode 
corresponding to the wave-vector q 2 occurs only when a coupling between two neighbouring 
chains is established through vertical contact. To identify the threshold for the transition 
from the modulated liquid to modulated crystal phase, we have plotted dpd/dcj) as a function 
of cp in Fig. 0(b). It can be seen that dpa/dcfi has a pronounced peak at certain cj> signifying 
sharp change in pd at that value of cj>. We identify the peak in dpd/dcj) as the onset of 
the freezing transition leading to a non-zero density mode pd- These area fractions for the 
zero-field and infinite-field are 0.702 and 0.688, respectively. 

In Fig. 0, we have plotted the variation of pd and dpd/dcj) with changing area fractions for 
some typical values of field strength (f3V = 0.1, 1.0, 10.0,2000.0, where (3 = l/fc^T). From 
the plots we have identified, for each value of (3V Q1 the corresponding value of cj> where the 
onset of freezing occurs. Using these values, we have thus constructed the phase diagaram 
of the hard-sphere colloidal system in the presence of an external laser field, which is plotted 
in Fig ||. 

As noted above, in the absence of the laser field, the onset of freezing takes place at an 
area fraction value of 0.702 . When the field is switched on, the density modes (pi) corre- 
sponding to the modulating wave- vector become nonzero, and this facilitates the deveopment 
of the the other density modes (pd) causing the freezing transition to happen at a smaller 
cj>. As we can see from the phase diagram, the value of cj> required for freezing continues to 
decrease till a (3V value of 2.0, where it is down to 0.667. Then it increases again to reach 
the infinite-field limit of about 0.69. A possible explanation for this latter increase could be 
the following. As f3V D gets increased to values much larger than unity, the motion of the 
particles become more and more restricted to be around the lines q D y = (2n + l)n. Thus, 
beyond a threshold value of the field strength, the influence of particles of neighbouring lines 
on the motion of particles on a particular line decreases. This manifests in the fact that the 
condensation of the density mode pd takes place at a higher area fraction compared to the 
values corresponding to field strengths below that threshold. With increasing field strength, 
the value of the critical cj> goes on increasing until it reaches the limiting value of 0.688. 

If we examine the phase diagram in Fig. |3|, we can see that for any fixed value of cj> 
less than 0.66, the colloidal system continues to be in modulated liquid state for all values 
of f3V Q . For values larger than 0.7, the system remains in crystalline form throughout the 
entire range of field strengths. The intermediate region is interesting. For example, consider 
a cj> value of 0.675. For low field strength, at this volume fraction, the system behaves like a 
modulated liquid. If one exceeds a certain field strength ( j3V Q = 1.4 ), the system undergoes 
a transition to a modulated crystalline state. But beyond j3V Q = 31, the hard sphere system 
again becomes a modulated liquid. Thus, from the measurement of the translational order 
parameter in Monte Carlo simulations, we find a possibility of re-entrant behaviour in the 
two-dimensional hard sphere system in the presence of an external modulation potential. 

To examine the three phases at cj> = 0.675, we have measured in our simulations the 
equilibrium real-space particle density < p(x, y) > for three different field strengths, as 
shown in Fig|4|(a)-(c). After equilibrating for 10 5 configurations, < p{x,y) > was measured 
by averaging over 1.5 x 10 5 configurations. From the plots, we can see that for (3V equal 
to 0.1 , the system is in a liquid state, the modulations not being very prominent. When 
(5V is equal to 2.0, the system is nearly crystalline and for j3V Q = 1000 , the system loses its 
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crystalline structure and again becomes a modulated liquid. This confirms the occurrence 
of a re-entrant melting as a function of the strength of the modulating potential. This is 
consistent with the recent prediction by Frey et al I|, but it remains to be shown that the 
melting here is a continuous transition and is mediated by the unbinding of dislocations. 
More extensive simulations are needed to demonstrate the latter. Meanwhile, it would be 
interesting to carry out experiments involving hard sphere colloidal particles and observe 
whether such a reentrant melting occurs as the strength of the laser field is increased. 

In conclusion, we have studied the effect of the external laser field modulation on a 
two-dimensional system of colloidal particles interacting via a hard sphere potential. By 
identifying the maxima in the rate of variation of the order-parameter ( as a function of 
the packing fraction ) as the onset of condensation of the density mode, we have constructed 
the phase diagram for the system. For a certain range of packing fractions, re-entrant melting 
occurs in the system as the strength of the field increases. Measurement of real-space density 
for (p = 0.675 confirms that the system first freezes into crystalline form but later on becomes 
a modulated liquid on increase of the field strength. A method has been proposed recently 
to calculate isothermal elastic constants for the hard disk triangular solids in two dimensions 



using Monte Carlo simulations |TT |. It would be of interest to check (using this method) 
whether these elastic constants soften when the field strength is increased beyond a certain 
value in our systems. 

After our work was completed, we learnt of an independent and more extensive Monte 
Carlo simulational study on the same system by Strepp et al jl2). Using the the cumulant 



intersection method, they have obtained a phase diagram has been obtained which is quali- 
tatively and quantitatively similar to the one obtained by us. This further strengthens the 
case for the existence of the re-entrant modulated liquid phase in a hard-sphere colloidal 
system subject to an external laser field modulation. 
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FIGURES 

FIG. 1. Variation of pd and dpd/dcj) for (3V = and for infinte field . 

FIG. 2. Plots of pd and dpd/dcf) for some typical values of field strength, 

(3V = 0.1,1.0,10.0,2000.0 . 



FIG. 3. Phase diagram of the colloidal system for 900 particles in an external laser field. 

FIG. 4. Contour plots of average density for cf> = 0.675 in presence of external field, (a), (b) 
and (c) correspond to (3V Q = 0.1,2.0, 1000.0 respectively. 
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Fig.4b 
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<p(x,y)> for p V q = 1000 
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Fig. 4c 
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